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Abstract. We have proposed a model to analyze the growth kinetics of lysozyme 
crystals/aggregates under non- isothermal conditions. The model was formulated 
through an analysis of the entropy production of the growth process which was obtained 
by taking into account the explicit dependence of the free energy on the temperature. 
We found that the growth process is coupled with temperature variations resulting in 
a novel Soret-type effect. We identified the surface entropy of the crystal/aggregate as 
a decisive ingredient controlling the behavior of the average growth rate as a function 
of temperature. The behavior of the Gibbs free energy as a function of temperature 
is also analyzed. The agreement between theory and experiments is very good in the 
range of temperatures considered. 
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1. Introduction 

Experimental and theoretical studies on the active role of temperature in protein crystal 
formation and other related problems of soft-matter aggregation could open a relevant 
possibility toward controlling the morphological and structural properties of protein 
crystals through an adequate tunning of the thermal conditions in which they grow 
[U |2j [31 IH O Ej. The importance of this research becomes manifest in biomedical jl] 
and clinical studies where the effect of temperature on forming certain (disease-assisted) 
plagues or fibrils was observed to help in partitioning aggregates into several groups of 
pathomorphological interests [7]. 

The research on the active role of temperature in exploring optimal nucleation 
and growth conditions of protein crystals or aggregates is mostly empirical. A large 
number of experiments conducted in thermostatted conditions had contributed to the 
comprehension of the sensitive-to-temperature kinetics by monitoring the temperature 
dependence of related thermodynamic quantities such as the growth rate of the 
crystals/aggregates [H [21 El HI E|- For example, in recent years the effect of 
temperature on average (110) face growth rates in the range of 277 — 293K at 
different aqueous solution and salt concentrations has been measured carefully [Hj • This 
experimental study has shown that the crystal formation appears ultimately to be a 
common product of interfacial and solution involving mechanisms. The formation, seen 
in terms of interface dynamics, involves the existence of the so-called depletion zone [1] , 
that is, a screened interfacial force field of electrostatic and hydrophobic nature whose 
characteristic length is temperature dependent O El IH1 [TO] - 

Here, we present a theoretical model accounting for the average growth rate of 
crystals/aggregates as a function of temperature. We start by calculating the mesoscopic 
entropy production [11] of the growth process that contains contributions associated 
to two variables, the radius and the temperature of the crystal/aggregate. These 
contributions lead to a novel Soret-type effect coupling the dynamics of the radius with 
the variations of the temperature, in a way similar to that of the conventional thermal 
diffusion [TJl [131 EH HSj • After explicitly considering the dependence of the Gibbs free 
energy on temperature through the surface energy term, a formula is obtained and 
used to make quantitative comparisons with experiments. By discriminating clearly the 
effects of the surface energy and temperature, our model may also help in discerning the 
degree of importance of the electrostatic screening effect and of temperature variations 
on the appearance of different surface morphologies and therefore of the output forms 

pq hsi nzi na m ess eu. 

The proposed model may constitute a comprehensive and useful theoretical 
framework for describing the formation of lysozyme crystals/aggregates and its growth 
in temperature-controlled conditions by addressing: a) The free-energy behavior as 
a function of temperature in lysozyme crystallization, and b) how this behavior is 
influenced by the temperature through the contribution of the surface entropy of the 
crystals/aggregates. Finally, c) we formulate a quantitative model accounting for the 
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behavior of the crystal growth rate as a function of temperature that takes into account 
a coupling term associated to the thermal expansion of the system and mediated by the 
Soret-type effect. The results of the model show good agreement with experiments. 

The paper is structured as follows. In Sec. 2 we calculate the mesoscopic entropy 
production leading to the Soret-type diffusion in cluster size _R-space. Then, in Sec. 3 
we discuss the dependence of the Gibbs free energy on temperature and its influence in 
crystal's surface energy and entropy, and show how the Soret-type coefficient and the 
average growth rate of crystals/aggregates emerge from our analysis. Sec. 4 is devoted 
to formulate the model for the average growth rate and to compare with experiments. 
In the final section, we shall conclude on the main findings and their inherent effect on 
the model soft-matter aggregating system. 



2. Entropy production in protein crystal formations: Size and temperature 

The characterization of the nonequilibrium state of the protein aggregate is performed 
through the knowledge of the probability density P = P(R,t), of finding the system 
with an effective radius R G [R, R + dR] at time t. This quantity obeys the continuity 
equation 

^P(R,t) + M J(R,t) = 0, (1) 

where J(R, t) is an unknown probability current in i?-space. Protein aggregation is an 
irreversible process and as such is governed by the thermodynamics laws. To determine 
the value of J(R,t), we will make use of the Gibbs' entropy postulate in the form [11] 

S(t) = S leq -k B [pin-^dR, (2) 
J ±\ eq 

where Sit) and S\ eq are the nonequilibrium and local equilibrium entropies of the system. 
The equilibrium value of the probability Pi eq (R) is given by the canonical distribution 
function PU [221 1 



Pleq(R) = ^e- W ™' kBT . (3) 

Here, Z = Z(T) is the canonical partition function that depends on the temperature of 
the aggregate T which is imposed by a heat bath. W m i n is the minimum reversible work 
required to create the actual state of the aggregate and ks is the Boltzmann constant. 
It is important to notice here that the radius R of the aggregate is in general a function 
of the temperature, R = R{T). This fact will be taken into account in what follows. 

At constant pressure and temperature, the minimum-work associated to the 
formation of the protein aggregate can be identified with a free energy difference [11] 

W min (R,T) = AG(R,T). (4) 

Protein aggregation entails rearrangements of its intrinsic structure, a process which 
is only possible at a certain energetic cost. The rate at which energy is dissipated is 
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To~, where a is the entropy production rate [TT] . This quantity follows from the time 
derivative of Eq. ([2]) after taking into account that in the local equilibrium state 



dS leq (t) _ J ,_ 1 dE{t) 



dJ n 



- T - l r^dR (5) 
dt dt J OR y 1 

where J q (R,t) is the amount of heat exchanged per unit time during the growth of the 

crystal/aggregate [24] . An integration by parts after employing (pQ) yields an evolution 

equation for the local entropy s(R, t) of the form: ds/dt + 3J s /dR = a. Here, J s is the 

entropy current and the entropy production rate a(R,t) is given by 



(6) 



J q dT J\{ 8AG \ k B T dP 
a{ ' } ~ ~T^dR ~ f [\dRJ T + 

This equation establishes that energy dissipation during the evolution in time of the 
system is due to two contributions. The first one comes from the dependence of the 
radius of the aggregate on the temperature. The second contribution comes from 
the driving force leading the aggregation in i?-space, and coming from the explicit 
dependence of the nonequilibrium chemical potential /i(-R, T, t) on the radius of the 
crystal/aggregate. Here, fi(R,T,t) is defined by 

fi(R, T, t) = fx leg (T) + k B T In [P/Pi eq ] , (7) 

where ni eq {T) = — k B T\n Z(T). In deriving Eq. §6§ we have introduced the modified 
heat current J q = J q — hfj that incorporates the enthalpy of formation of the 
aggregate hf and we have separated the variation of the chemical potential in the form: 
dfi = (dfi/dR) T p dR+ (dfi/dP) T R dP + (dfi/dT) R p dT. In this relation the subindexes 
indicate the variables that should be considered constant when performing the derivative 

From the entropy production ([6]) it is possible to formulate the following linear law 
for J(R,t): 



J(R,t) 



"7 



P 



( dAG \ 
\~dRJ 



+ k B T 



T 



dP 
dR 



-D T P 



.dT 
dR' 



(8) 



where 7 and Dt are, respectively, the mobility and thermal diffusion-type coefficients in 
the space of crystal/aggregate radius R. After substituting Eq. ([8]) into ([T|) we obtain 
a Fokker-Planck type equation 



d_ 

Wt 



P(R,t) 



d_ 
dR 



7 



(dAG\ 
\dR~J, 



P + k B T~i 



dP 
dR 



+ 



d_ 
dR 



D T P 



dT\ 
dR 



(9) 



This equation constitutes the starting point in the formulation of our model accounting 
for the dependence on temperature of the growth rate of protein-type aggregations as, 
for instance, lysozyme aggregations in aqueous solutions [27J [281 129] . This dependence 
has been reported in experiments that confirm the sensitivity of the growth of the 
aggregates to temperature conditions [H |2l [31 IH [H |9] , and can be analyzed through an 
adequate choice of the Gibbs energy and the specific physical properties of each protein 
crystal/aggregate such as thermal expansion coefficient and surface energy. 
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3. Temperature dependence of the Gibbs energy and the Soret-type effect 

In order to quantify the effect of temperature in protein aggregations it is convenient to 
first derive the evolution equation for the average radius of the aggregate R(t), defined 
by the relation: R(t) = J RPdR. 

Thus, according to its definition, the expression for the average growth rate of the 
(spheroidal) aggregate V a = dR{t)/dt is obtained by multiplying Eq. (j9]) by R and 
then integrating by parts over i?-space. Assuming that fluctuations are small, then the 
probability can be approximated by a Dirac delta function of the form: P ~ 5(R(t)—R). 
One obtains 



dR(t) D \ 1 



dt I ksT 




(10) 



R 



where we have defined D = ksT^f and the Soret-type coefficient St = D T /D, in a way 
fairly reminiscent of that presented in [151 US] f° r conventional thermo diffusion. Eq. ([TO]) 
shows that under non-isothermal conditions, the growth rate of the aggregate depends 
on the diffusion constant and free energy barrier, as occurs in the isothermal case, and on 
local variations of the temperature along the radius of the crystal/aggregate. This new 
contribution results from the coupling between the growth kinetics with the temperature 
inhomogeneities that define a Soret-type effect in R-space. It is opportune to note that 
V a , as an averaged measure proposed to quantify the growth rate, may differ from the 
instantaneous growth rate proposed earlier [22j [18] for the aggregations exploited as 
isothermal formations (301 EJJ- 

At equilibrium (V a = 0) the Soret-type coefficient is determined from Eq. f[T0[) by 
the dependence of the free energy on the temperature in the form 

St = — ~; — ~ 

k B T 

This relation is similar to that reported in the literature for the Soret effect examined 
in particles' position space [131 E] • 

In this way, the explicit dependence of St on the temperature can be obtained by 
assuming that the free energy change is given by the expression of the classical nucleation 
theory 

AG = -(4<K/3v p )R 3 Afl + AttTR 2 , (12) 

where Ap,(T) = jl c — is the electrochemical supersaturation given by the 
electrochemical potential difference between the crystal jl c and the protein solution 
/ij. Here, T is the surface energy. This expression assumes the presence of a collective 
interfacial tension mechanism for the aggregation, originally proposed in Ref. [32] for 
some inhomogeneous protein (lysozyme) solutions. Moreover, we have tacitly supposed 
that in first approximation the form of the protein can be assimilated to an effective 
sphere with a gyration radius. Recent results on computer simulations of the structure of 
certain proteins under non-isothermal conditions [33] have shown that the fluctuations 
of their radius of gyration are small thus indicating the validity of such an effective 
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model. The symbol v p thus stands for the protein effective volume and the number of 
proteins in the aggregate can be obtained through the relationship n a = 4nR 3 /3v p . 

The temperature dependence of the physical parameters entering Eq. ( fl2|) can 
be inferred on the following basis. According to the classical nucleation theory, the 
chemical potential difference in Eq. ( Ti~2l is customarily assumed to be of the form: 
A/i(T) = kT In \C/ci\, with q the density of proteins in solution and C the density of the 
crystal/aggregate. In this approximation, m|C/cj| can be considered as independent of 
R and T within the range of temperatures of interest for the (dis)orderly aggregation 
such as the one of lysozyme type in aqueous solution, taking mostly place at pH ^4.5 
and up to about T ~ 310K (high-temperature aggregation) [191 E2]- In addition, one 
has to consider the dependence of the protein volume v p on the temperature. According 
to Ref. [34] . this dependence is negligible for proteins in their native conformations and 
relatively weak for denatured proteins over the range of temperatures of interest. For 
convenience, in what follows we will assume that the protein volume is constant. 

As a consequence of the previous considerations, we may assume that the only 
parameter depending explicitly on the temperature in Eq. f fl2|) is the surface energy 
r(T). It may be expressed as a power expansion of the temperature with respect to a 
characteristic temperature T . Up to second order one obtains 



where T is the value of the surface energy at T = T . This relation can be written in 
terms of the surface entropy [35]: S s = — (dT/dT) p . Substituting the previous relation 
into ( fl3l) we finally obtain 



where b = (l/2)(d 2 T/dT 2 ) p . Here, the characteristic value of the surface entropy at 

(s) 

T = T is assumed to be positive: Sq > 0. Concerning the value of the parameter 
b two possibilities can be analyzed [25]: i) For b < 0, the free energy is a monotonic 
decreasing function of T. ii) For b > 0, after a first linear decrease which is dominated 
by the volumetric and surface entropy terms, the free energy presents a minimum due 
to the quadratic contribution. This second possibility reproduces the expected behavior 
of the free energy as a function of the temperature [34J, as shown in Fig. 1 obtained 
with the use of Eqs. (114]) and (11 2p . The appearance of this minimum makes it possible 
the control of the growth dynamics, as will be shown in the analysis of the behavior of 
the growth rate. Figure 2 shows in turn the Soret-type coefficient (fTTj) as a function of 
the temperature. It is interesting to mention that this coefficient exhibits a change of 
sign similar to that observed in the conventional Soret effect [121 EH EEH US] . The four 
cases shown in Figs. 1 and 2 correspond to the values of the parameters obtained from 
the fit shown in Fig. [3] and given in Table [1] and in the text. The value of To in each 
case corresponds to the initial point of the growth rate in Fig. |3j 

Here, it is worth mentioning that for the long stage isothermal growth of the 
crystal/ aggregate, it was proposed in Refs. pUl E21 EQl [31] that the free energy AG(R, T) 




(13) 



r(T) ~ r - S { S \T - T ) + b(T - To) 2 , 



(14) 




Figure 1. The free energy given by Eq. (|T!21) as a function of the temperature in 
the case when the surface energy T(T) is given by Eq. (TH| . The free energy has the 
expected behavior, cf. Ref. [33]. The values of the parameters used are those obtained 
in the fit shown in Fig. [3J 



can be modeled through the near-interface electrochemical supersat uration of the system 
by AG(R, T) = k B T In \5i\, where 5i = 5i(R,T) = Ci(R,T)/[C — Ci(R,T)] is the inverse 
near-interface solubility Tii(R,T) = 1/Si. This quantity is expected to be measurable 
in terms of the inverse of the corresponding depletion-zone supersaturation. Here C is 
again the density of the aggregate and now Ci(R, T) is the near interface concentration of 
proteins, whose expression is given in turn by Ci(R,T) = c [l + {2Y / R){1 / ksT)]. Here, 
Co is the concentration of proteins in the solution upon flatness of the aggregate surface 
and Ci(R,T) is viewed as the near interface density that incorporates a dependence on 
the curvature and temperature. Notice that now T/ksT becomes the Gibbs-Thomson 
capillary length [22]. The dependence of Ci(R,T) on T introduces a correction of the 
supersaturation making it of non-ideal nature. 

This protein aggregation model makes it possible to incorporate second and higher 
order contributions of a virial expansion of the local chemical potential. The dependence 
of the second virial coefficient on temperature has been analyzed in Ref. [10]. Finally, 
it is worth to mention that the generalized form of the near interface concentration 
Ci(R,T) is consistent with the Gibbs-Thomson boundary condition coming from Eq. 
(I12p when AG is minimized. 

4. Comparison with experiments 

From the previous analysis, one may conclude that in non-isothermal experiments 
performed at constant supersaturations [21 [8], the Soret-type effect characterized by 
the coefficient St can be relevant in determining the dependence of the growth rate V a 
on the temperature. 
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Figure 2. The Sorct-typc coefficient St given by Eq. (fTTj) combined with p^|) and 
(|14[) as a function of the temperature and corresponding to the free energy of Figure 
1. The values of the parameters were those obtained from the fit in Fig. [3] 



The behavior of the average growth rate for lysozyme aggregation in terms of 
temperature at constant supersaturation was measured in Refs. [2] and [8]. The 
experimental data (symbols) are shown in Fig. [3] together with the results (lines) given 
by the expression 

D 



Va 



^AG(R(t),T) 



J T 



I^H^l- (15) 

obtained from our model after substitutiing Eq. ( TTT1) into ( flOl) and invoking the 
definition of the thermal expansion coefficient oct{T) = (1/V) [dV/dT] , which in the 
present case leads to the relation «t = (3/R) [dR/dT] , [35] . 

This model relation for the average growth rate involves the competition of 
two forces coming from the dependence of the Gibbs energy on the radius and the 
temperature. The first term is the usual isothermal contribution whereas the second 
term comes from the temperature dependence of the Gibbs energy and is therefore 
influenced by variations the surface entropy S s and the thermal expansion coefficient 
olt- The results obtained lead us to conclude that a consistent description of the process 
necessarily entails to consider the coupling of the growth kinetics with the temperature 
through the second term in Eq. (fT5]) . 

In order to reproduce the experimental data, we have assumed that near the 
characteristic temperature To the thermal expansion coefficient can be expanded in 
the form olt — «o + $ [{T — T ) /T ] a , with a reference value of the thermal 
expansion coefficient. The scaling prefactor 5 is related to the derivative of the thermal 
expansion coefficient with respect to the temperature and a is a scaling exponent. 
From this expression for ar one can therefore infer that the dependence of the average 
crystal/aggregate radius on temperature is of the form R(T) = Rq exp[J-"(T)] with Rq 
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Figure 3. Average growth rate in terms of temperature (°K) for four different 
values of the supersaturation corresponding to protein concentrations of 35mg/ml 
(solid line), 30mg/ml (dashed line), 20mg/ml (dash-dotted line) and 10mg/ml (dash 
double-dotted line). The symbols are experimental results taken from [5] whereas the 
lines correspond to our proposed fitting model, Eq. (fTrJ)) . see also (|15p . The dotted 
line was obtained using the same values of the parameters as the dash-dotted line 
but neglecting the contribution of the second term at the right hand side of (|15|) . 
The difference emphasizes the fact that the new effect may be important for correctly 
accounting for the data at large temperatures. The values of the parameters are given 
in Table [T] and in the text. The value of the surface energy for the case of the dash 
double-dotted line was Tq — 4.5gs~ 2 . 



the value of the radius at To and T{T) = ff ctrdT. For simplicity in the comparison 
with experiments [8] , we have assumed the following linearized expression for the radius: 



R{T) ~ #0 [1 + "o( T - T o) + 5T 1+e (T - T )- £ J where e = 1 + a. Finally, assuming that 
in the range of temperatures considered the phenomenological coefficient D is constant, 
we obtain the formula 



D 



k B T 



r + AT ( bAT 



2b AT - S, 



.^ii 



1 

Taj 



a T (T) 
C 



R{T) 



kT In 



R{Tf 



(16) 



where we have defined AT = T — T . 

The values of the parameters considered in the cases studied are the following: 
Rq = 0.1/xm, Up = 0.049 fim,5 = 0.123 and e = 0.37. The comparison was done 
using the value of the thermal expansion coefficient of lysozymes reported in Ref. [36J : 
ao = 1.35 ■ lO -4 ^ -1 , valid whithin the temperature range: 298 — 323 K. The value 
obtained for the surface energy: To = 5.975gs -2 is compatible with that of the surface 
tension in a lysozyme solution reported in Ref. [3?]: 7 ~ 35gs~ 2 . For the diffusion 
coefficient we have obtained D ~ 1 • 10~ 10 /^m 2 s _1 . This apparently small value may 
be understood by considering that the diffusion coefficient in solids can be expressed as 
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Table 1. Values of the parameters entering Eq. (|T6| as obtained by fitting the 
experimental data shown of Fig. [3J which are in aggreement with those given in Ref. 
E] and [35]. 

[35] : D = D exp(—E ac /kBT). Here D is a reference value of the diffusivity and the 
so called exponential pre-factor that depends on the characteristic energy E c . Taking 
into account that Dq ~ 2 • 10 -2 fim 2 s _1 , see Ref. [38], and assuming that T ~ 280K, 
then we may estimate that E c ~ 19A;bT. This result is compatible with the value of the 
adsorption energy reported in Ref. [38], E ads ~ 12/c B T (see also Ref. [2]). Finally, the 
values of the remaining parameters are given in Table [TJ 

From the results obtained, we can conclude that the maximum of the average 
growth rate V a (T) is due to the presence of a minimum of the Gibbs energy and 
therefore a consequence of the competition between volumetric and surface effects. For 
a fixed value of D, the initial steep and the amplitude of the maximum of V a are 
associated respectively to the supersaturation parameter and to the surface entropy of 
the crystal/aggregate, that dominates the surface energy at low temperatures, this fact 
is shown in the first and second columns of Table [TJ The second term of the surface 
energy controls in turn the decrease of V a for larger values of the temperature. For 
smaller values of the parameter b (having the dimensions of a heat capacity per unit 
area), the tail of the curve becomes larger thus improving the fit, see Fig. 3 and Table 
[TJ The decay of V a is also affected by the value and dependence on temperature of 
the thermal expansion coefficient otT- This fact is illustrated in our main finding, Fig. 
3, through the dash-dotted and dotted lines. The dash-dotted line was obtained by 
considering the non-isothermal correction from expression (1151) or ([TBI , whereas the 
dotted line is the result in the isothermal case when the free energy of Eq. (fT5l) does not 
depends on temperature. It is convenient to emphasize that the contribution coming 
from the Soret-type effect may be relevant in explaining experimental results as those 
discussed here, because the errors reported for the measurements of V a (not shown here) 
are significantly larger for intermediate values of T (i.e., around the maximum) whereas 
they were small for the high temperature regime, see Ref. [8]. A similar behavior of the 
average growth rate in terms of the temperature has been also reported from simulations 
in Ref. [39] • 

At constant supersaturation, the results for the ratio 47^/1/31^0 shown in the 
third column of Table [TJ indicate that when the surface energy Tq of the crystal/aggregate 
is smaller (solid and dashed lines) the maximum of V a (T) appears at temperatures larger 




Figure 4. Standard crystallization enthalpy AH(T) as a function of temperature 
coming from Eqs. (fT2|) and (|T7|) for the values of the parameters obtained from the 
fitting shown in Fig. [3J 



than those corresponding to the case of smaller surface energies, in which the maximum 
decreases in amplitude and shifts towards the region of lower temperatures (dash-dotted 
and dotted lines). Surface energy inhibits the formation of the crystal/aggregate whereas 
surface entropy favours it, see the fourth column in Table [TJ 

Another relevant and measurable quantity [TOl [9] that can be obtained from our 
model is the enthalpy change with temperature. The standard crystallization enthalpy 
change denoted by AH can be determined from the Gibbs-Helmholtz equation. The 
sign of AH indicating endothermic (AH > 0) or exothermic (AH < 0) crystallization 
can be obtained from [9j [10] 

Making use of Eqs. f fl2|) . f JT^|) and ffTTl) with the values of the parameters obtained in 
the fit of the experimental data, we obtain the behavior of the enthalpy in terms of 
temperature shown in Fig. HI As in the case of the Soret-type coefficient, the change of 
sign of the enthalpy coincides with the presence of a minimum of the free energy shown 
in Fig. 1. The quantity AH was calculated in [10] for the case of oxidized-and- mutant 
hemoglobin C crystal formation. 



5. Conclusions 



We have analyzed the influence of temperature variations on protein crystallization, a 
phenomenon significant for example in biomedicine where small changes in temperature 
from one tissue to another or in the same tissue with time may have an important 
influence on the crystallization of proteins that occur in living systems [U,L7]. 



Controlling protein crystal growth rate by means of temperature 



12 



In particular, we have presented a kinetic model which considers the dependence of 
the Gibbs free energy on temperature to calculate the average growth rate of lysozyme 
crystals. An analysis of the mesoscopic entropy production of the growth process 
reveals the presence of two contributions to the growth rate. The first one is similar 
to that for isothermal growth and can be computed from classical nucleation theory; 
the second stemming from the coupling between growth, viewed at short time scales 
as diffusion along the size of the cluster coordinate, and temperature differences. This 
second contribution constitutes a novel Soret-type effect. As we have shown, this effect 
is essential to correctly explain experimental results for the average growth rate of 
crystals/aggreagates in terms of the temperature. 

We used our model to reproduce experiments performed at constant supersaturation 
in the high temperature range. The quantitative agreement between experiments and 
our model formula ( Tl6|) is very good even at higher temperatures, that are difficult to 
capture with a single expression for the free energy. Other models based on the kinetic 
roughening hypothesis also reproduce the observed crossover of the growth rate (2J HO]. 
The obtained fit for this quantity is good at low temperatures but at high temperatures 
it deviates substantially from the measurements. In this domain, the fit provided by 
our model is better. The fact of reproducing the growth rate in the whole temperature 
range allows for a systematic interpretation of the experiments in terms of the surface 
energy and entropy of the crystal/aggregate. Our results also show the important role 
played by the thermal expansion coefficient in determining the behavior of the average 
growth rate with temperature. 

The model can be generalized to describe non-spherical, polygonal-type crystalline 
nuclei by explicitly considering the local curvatures R~ l of their z-faces. In this case, each 
face will show different behaviors of the corresponding surface energies To,*; entropies 
S St i and components of thermal expansion tensor leading to different behaviors of 
the corresponding growth rates V a ,i(T), [H El El EH H2]- 

The present study may contribute to a better understanding of several experimental 
efforts [21 [3J EJ [10] examining the dependence of pure crystal/aggregate morphologies 
with temperature. Future examination of the growth process involving nucleation 
promoters such as salts is left for another study [H [29]. 
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